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ABSTRACT 



Context. Observing supernova remnants (SNRs) and modelling the shocks they are associated with is the best way to quantify the 

energy SNRs re-distribute back into the Interstellar Medium (ISM). 

Aims. We present comparisons of shock models with CO observations in the F knot of the W28 supernova remnant. These comparisons 

constitute a valuable tool to constrain both the shock characteristics and pre-shock conditions. 

Methods. New CO observations from the shocked regions with the APEX and SOFIA telescopes are presented and combined. The 

integrated intensities are compared to the outputs of a grid of models, which were combined from an MHD shock code that calculates 

the dynamical and chemical structure of these regions, and a radiative transfer module based on the 'large velocity gradient' (LVG) 

approximation. 

Results. We base our modelling method on the higher J CO transitions, which unambiguously trace the passage of a shock wave. 

We provide fits for the blue- and red-lobe components of the observed shocks. We find that only stationary, C-type shock models can 

reproduce the observed levels of CO emission. Our best models are found for a pre-shock density of lO'' cm"', with the magnetic field 

strength varying between 45 and 100 yuG, and a higher shock velocity for the so-called blue shock (~25 km s"') than for the red one 

(~20 km s"'). Our models also satisfactorily account for the pure rotational H? emission that is observed with Spitzer. 



Key words. ISM: supernova remnants - ISM: individual objects: W28 - ISM: kinematics and dynamics ■ 
cesses: shock waves - Submillimeter: ISM - Infrared: ISM 



Physical data and pro- 



1. Introduction 

The interstellar medium (ISM) is in constant evolution, ruled 
by the energetic feedback from the cosmic cycle of star forma- 
tion and stellar death. At the younger stages of star formation 
(bipolar outflows), and after the death of massive stars (SNRs), 
shock waves originating from the star interact with the ambi- 
ent medium. They constitute an important mechanical energy 
■ input, and lead to the dispersion of molecular clouds and to the 
[ compression of cores, possibly triggering further star formation. 
. Studying the signature of these interactions in the far-infrared 
• and sub-mm range is paramount for understanding the physical 
and chemical conditions of the shocked regions and the large- 
scale roles of these feedback mechanisms. 

Supernovae send shock waves through the ISM, where they 
successively carve out large hot and ionised cavities. They subse- 
quently emit strong line radiations (optical/UV), and eventually 
interact with molecular clouds, driving lower-velocity shocks. 
Similar to their bipolar outflow equivalents, these shocks heat, 
compress, and accelerate the ambient medium before cooling 
down through raolecu lar emission (van Dishoeck et al. 1993, 
lYuan & NeufeIdll20TTl hereafter Yl 1). 

V aluable information has be en provided by 

ISO (ICesarskvetaLL 119991: ISneU et al.L l2005h and Spitzer 
dNeufeld et al.1 l2007l hereafter N07), but neither of those 
instruments provided sufficient spectral resolution to allow for a 



detailed study of the shock mechanisms. High-7 CO emission is 
one of the most interesting diagnostics of SNRs. CO is indeed 
a stable and abundant molecule, and an important contributor to 
the cooling of these regions, whose high-frequency emission is 
expected to be a 'pure' shock tracer Observations of the latter 
must be canied out from above the Earth's atmosphere. As part 
of a multi-wavelength study of MHD shocks that also includes 
Herschel data, we present here the first velocity-resolved CO 
(11-10) observations towards a prominent SNR-driven shock 
with the GREAT spectrometer onboard SOFIA, and combine 
them with new lower-7 ones in a shock-model analysis. 



2. The supernova remnant W28 



W28 is an old (>104^ vr. lClaussen eTaPl 19991) SNR in its radia- 
tive phase of evolution, with a non-thermal radio shell centrally 
filled with thermal X-ray emission. Lying in a c omplex region of 
the Galactic disk at a distance of 1.9+0.3 kpc (Velazque z et al.L 
2002), its structure in the 327 MHz rad io continuum rep resents 
a bubble-like shape of about 40x30 pc (iFrail et al.Lll993h . Early 
on, molecular line emission peaks, not associated with star for- 
mation activity, but revealing broad lines, were suggested as ev- 
idence for interaction of the remnant with surrounding molecu- 
lar clouds ( Wootten , 1981 ). Later studies spati ally resolved the 
shocked CO gas layers from the ambient gas ( Frail & MitcheUi 
119981: lArikawa et"aLL Il999h . OH maser spots line up with flie 
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post-shock gas layer s (iFrail et al.l Il994t IClaussen et all Il997t 
iHoffman et all 120051) . for which the strongest masers VLBA po- 
larisation stu dies yield line-of-sigh t magnetic field strengths of 
up to 2 mG dClaussen et all Il999h . Pure rotational tran sitions 
of H2 have been detected with ISO dReach & Rhd I2OOQ) and 
were more recently observed with Spitzer, better resolved spa- 
tially and spectrally, by N07 and Yll. 

Recently, very high energy (TeV) y-ra y emissi on has been 
detected by HESS dAharonian et all l2008l) . Fermi (lAbdo et all 
l20ltf ). and AGILE dGiuliani et all 1201 Ol) . spatiallv slightly ex- 
tended and coincident with the bright interaction zones, W28-E 
and -F. If interpreted as the result of hadronic cosmic ray in- 
teractions in the dense gas (7r° decay), a cosmic ray density en- 
hancement by an order of magnitude is required (which is sup- 
plied/accelerated by the SNR). 

2.1. Sub-mm CO observations of W28F 



APEXQ dGiisten et all |2006|) observations towards W28F were 
conducted in 2009 and will be the subject of a forthcoming pub- 
lication (Gusdorf et al., in prep.). For the present study, we used 
100"xlOO" maps in the '^CO (3-2), CO (3-2), (4-3), (6-5), and 
(7-6) transitions, described in Appendix lAl 

In Fig. [U the velocity-integrated CO (6-5) broad-line emis- 
sion of W28F is shown overlaid with the CO (3-2) emission 
(white contours): a north-south elongated structure of about 
100" height and 30" width traces the same warm accelerated 
post-shocked gas. In our high-resolution CO (6-5) data the struc- 
ture is resolved, though probably still sub-structured similar to 
what is seen in H2, e.g., Yll. Comparison with the distributions 
of excited H2 and OH masers (whose locations also mark the 
leading edge of the non-thermal radio shell) suggests a textbook 
morphology of an SNR-molecular cloud interaction: the shock 
propagates E-NE into the ambient cloud that extends east for 
several arcmins. Hot H2 and OH masers mark the first signposts 
of the shock-compressed gas. Farther downstream, the gas cool- 
ing is seen prominently in warm CO. The shock impact appears 
edge-on, but the fact that high - projected - streaming veloci- 
ties are indeed observed (-30 km s ' with respect to the ambient 
cloud) requires a significant inclination angle. 

We have selected the most prominent position in the south- 
ern extension of W28F, marked with the blue circle in Fig. 1, for 
our shock analysis. This position was also covered by Spitzer, of- 
fering a set of complementary H2 data. Fig. |2] shows the APEX 
spectra obtained towards this position, in main beam tempera- 
ture units, all convolved to the 23.7" beam of the SOFIA ob- 
servations (Sect. I2.2I 1. Absorption notwithstanding, the spectra 
show Gaussian profiles, with all line wings extending in blue (- 
30 km s"') and red (-1-15 km s"') velocities with respect to the 
ambient cloud (~10 km s '). The higher the line frequency, the 
narrower the line profile, with a typical linewidth of about 20 km 
s ' for CO (7-6). We detect red-shifted line-of-sight absorption 
features in all lines up to CO (6-5), the deepest arising at the am- 
bient cloud velocity. OfF-position contamination results in minor 
absorption features at 20 and 25 km s"' in the (3-2) and (4-3) 
profiles. Comparison with the '^CO(3-2) profile is also shown in 
Fig. |2] An analysis of the line temperature ratio of '-CO/'^CO 
yields an optical thickness value of 3-7 in the wings of the '^CO 



100 




' This publication is partly based on data acquired with the Atacama 
Pathfinder Experiment (APEX). APEX is a collaboration between 
the Max-Planck-Institut fiir Radioastronomie, the European Southern 
Observatory, and the Onsala Space Observatory. 



offset (arcsec) 

Fig. 1. Overlay of the velocity-integrated CO (6-5) (colour back- 
ground) with the CO (3-2) (white contours) emission observed 
with the APEX telescope. For both lines, the intensity was in- 
tegrated between -30 and 40 km s"'. The wedge unit is K 
km s ' in antenna temperature. The CO (3-2) contours are 
from 30 to 160 cr, in steps of lOcr = 16 K km s"'. The half- 
maximum contours of the CO (3-2) and (6-5) maps are indi- 
cated in red and black, respectively. The dark blue circle in- 
dicates the position and beam size of the SOFIA/GREAT ob- 
servations. The APEX beam sizes of our CO (3-2), (4-3), (6- 
5), and (7-6) observations are also provided (upper right comer 
light green circles, see also Table lA.lb . The maps are centred 
at (R.A.[j2ooo]=18''01"'52^3, DeC[j2ooo]=-23°19'25"). The black 
and light blue hexagons mark th e po sition of the OH mase rs ob- 
served bv lClaussenetal.1 dl997h and lHoffman et all d2005h . 



(3-2 ), assuming a typical inte rstellar abundance ratio of 50-60 
(e.g.. lLanger & Penziasl[l993h . 



2.2. Far-infrared CO spectroscopy witti GREAT/SOFIA 

The observations towards W2 8F were con ducted with the 
GREA10 spectrometer (Hevminc Fet al.ll2012h during SOFIA'S 
flight from Stuttgart to Washington on September 21 2011. Only 
one position could be observed, towards the southern tip of the 
shocked cloud at offset (+7",-26") (Fig. [B- The CO (11-10) 
line was tuned to the frequency 1267.015 GHz LSB. The re- 
ceiver was connected to a digital EFT spectrometer (Klei n et al] 
2012) providing a bandwidth of 1.5 GHz with a spectral resolu- 
tion of 0.05 km s '. The observations were performed in dou- 
ble beam-switching mode, with an amplitude of 80" (or a throw 
of 160") at the position angle of 135° (NE-SW), and a phase 
time of 0.5 sec. The nominal focus position was updated reg- 
ularly against temperature drifts of the telescope structure. The 
pointing was established with the optical guide cameras to an 
accuracy of ~5". The beam width 0„ib is 23.7"; the main beam 
and forward efliciencies are //mb = 0.54 and rif - 0.95. The inte- 
gration time was 13 min ON source, for a final rm.s of 0.66 K. 



' GREAT is a development by the MPI fiir Radioastronomie and 
the KOSMA/Universitat zu Koln, in cooperation with the MPI fiir 
Sonnensystemforschung and the DLR Institut fiir Planetenforschung. 
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Fig. 2. CO transitions observed in the position (+7", -26") indi- 
cated in Fig.[T] APEX (3-2), black (corresponding '^CO, green); 
(4-3), pink; (6-5), dark blue; (7-6), Hght blue; and SOFIA (11- 
10), red. The '^CO (3-2) and CO (11-10) profiles were multi- 
plied by three for comparison purposes. Respective CO spectral 
resolutions are 0.212, 0.318, 0.159, 0.272, and 0.693 km s"', and 
0.664 km s' for '^CO (3-2). 

The data were ca Ubrated with the KOSMA/GREAT calibrator 
(iGuan et al.ll2012h . removing residual telluric lines, and subse- 
quently processed with the CLASS softwarqfl 

The CO (11-10) spectrum, overlaid on the sub-mm lines in 
Fig. |2] reveals a markedly different profile: weak emission only 
is seen from the high-velocity gas, while the line is more promi- 
nent at low velocities. The profile basically follows the shape of 
the optically thin '^CO(3-2). 

3. Discussion 

3.1. The observations 

Although most likely part of a single original shock clump, we 
separated the profiles into a blue lobe (-30 to 7.5-12.5 km s"') 
and a red lobe (7.5-12.5 to 40 km s '), and fited the data inde- 
pendently. The uncertainty of the upper (lower) limit of those 
ranges reflects our lack of knowledge of the ambient velocity 
component. In Fig. [3] we plot the velocity-integrated intensities 
of the CO lines, all convolved to the same angular resolution 
of 23.7", against the rotational quantum number of their upper 
level in a so-called 'spectral line energy distribution' (or SLED, 
filled black squares with errorbars). The underlying assumption 
is that the filling factor is the same for all CO lines, which is val- 
idated by the similarity between their emitting regions (see for 
instance Fig. [1] where the half-maximum contours of CO (3-2) 
and (6-5) coincide at the available resolution). For Fig. |3] the 
assumption was made that all observed transitions present a cir- 
cular emission region of radius 25", corresponding to a filling 
factor of 0.53 (compatible with our maps, see Fig. [1}. 

The upper (lower) panel shows the diagram associated with 
the blue (red) shock component. The lower limit to the integrated 
intensity corTesponds to the integration of the observed profile in 
the smallest velocity range, whereas upper limits were obtained 
by integrating Gaussian fits to the observation on the largest ve- 
locity range. The fits were adjusted to recover the shock flux lost 
through absorption and based on the un-absorbed parts of the 
profiles. Although yielding high erTorbars on our measurements, 
specially for low /up values, this method also provides the most 
conservative approach to our blue-red decomposition of the CO 
emission. 



http://www.iram.fr/IRAMFR/GILDAS 



3.2. The models 

We then compared the resulting integrated intensity diagrams 
to modelled ones. To build those, we used a radiative trans- 
fer module based on the 'large velocity gradient' (LVG) ap- 
proximation to characterise the emission from the CO molecule 
over outputs generated by a state-of-the-art model that calculates 
the structure of one-dimension, stationary shock layers (or ap- 
proximations of non-stationary layers). This method has already 
been used and extensively introduced in Gusdorf etal. (2008!). 
Since then, the LVG module has been modified to incorporate 
the latest colli s ional rate coefficients of CO with H2 computed 
by lYang et aP (l2010h . but the shock model is the same. Based 
on a set of input parameters (pre-shock density, shock veloc- 
ity, type, and age, and magnetic field parameter value b such as 
B[yuG] - b X -\Jnii [cm^^]), it calculates the structure of a shock 
layer, providing dynamical (velocity, density, temperature), and 
chemical (fractional abundances of more than 125 species linked 
by over 1000 reactions) variables values at each point of the 
layer. The relevant outputs are then used by our LVG module, 
which calculates the level population and line emissivities of 
the considered molecule. In the present case, the line tempera- 
tures were computed for transitions of CO up to (40-39) at each 
point of the shock layer, over which they were then integrated 
to form our model's SLED. Because of the importance of H2 
(key gas coolant, abundant molecule), its radiative transfer was 
treated within the shock model. Its abundance was calculated 
at each point of the shock layer, based on the processes listed 
in iLe Bourlot et all (l2002h . and the populations of the 150 first 
levels were also calculated inside the shock code, their contribu- 
tion to the cooling being included in the dynamical calculations. 
It is hence also possible to associate an H2 excitation diagram to 
each model of the grid, see Sect. IB. 21 

Our model grid consists of a large sample of integrated in- 
tensity diagrams, obtained for stationary C- and J-type, and non- 
stationary CJ-type shocks. It covers four pre-shock densities 
(from lO-' to 10'' cm"-'), 0.45 < b < 2, over a range of veloc- 
ities spanning from 5 km s"' up to the maximum, 'critical' value 
above which a C-shock can no longer be maintained, which de- 
pends on the other parameters. We then independently compared 
the observed blue- and the red-shock component to the whole 
grid diagrams. We used a x^ routine that was set up to provide 
the best fits to our purest shock-tracing lines: CO (6-5), (7-6), 
and (11-10), to reduce any effect induced by ambient emission 
contamination on the lower-lying transitions. 

3.3. The results 

The upper (lower) panel of Fig. |3] shows our best-fitting mod- 
els for the blue (red) component. We found that only stationary, 
C-type shock models can reproduce the levels of observed CO 
(6-5), (7-6) and (11-10) integrated intensities. In both cases, 
the pre-shock density is lO'* cm"^, for which a stationary state 
is typically reached for a shock age of lO'* years, in agree- 
ment with the age o f the remnant of 3.5x10"* years quoted by 
iGiuliani et alJ (l2010l) . The modelled shock velocities, 25 and 
20 km s~\ respectively, are rough upper limits to the observed 
ones, an expected conclusion given the one-dimensional nature 
of our models. Additionally, our models provide respective con- 
straints to the magnetic field component perpendicular to the 
shock layers, with parameter b of I and 0.45. The difference 
between those two values can be explained by a projection ef- 
fect. In both cases, the direct post-shock density is of the order 
of 2-3xl0^cm"^, yielding expected post-shock values for this 
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Fig. 3. Best-model comparisons between CO observations and 
models for each observed shock: the blue- (upper panel) and the 
red- (lower panel) components. Observations are marked by the 
black squares, our individual best-shock models are in blue line 
and circles (^i) or diamonds {(2), the warm layer that we used 
to compensate the ambient emission affecting the (3-2) and (4- 
3) transitions in green line and triangles, and the sum of each 
of these components is represented by a red line and circles or 
diamonds. 



component of 0.05-0. 55 mG, compatible with those inferred by 
IClaussen et al. ( 1999) (2 mG for th e total post-shock magnetic 
field), and bv Hoff^man et alj ( l2005h (0.3-1.1 mG based on the 
OH maser measurements in the positions indicated in Fig. \T\lo- 
cated within the blue beam). Finally, our models produce typical 
OH column densities of a few lO'^cm"^. This complies with the 
list of requirements set out by iLockett et al.l (1 19991) for exciting 
the observed OH 1720 MHz masers, together with the C- nature 
of the considered shocks, as well as the temperature and density 
conditions they create. 

The respective CO column density found in each of the red 
and blue shock models is 11 and 5xlO'^cm~^. However, our 
shock models failed to account for the important observed levels 
of CO (3-2) and (4-3) emission. We found it possible to model 
this lower-7 surplus emission by adding a thin layer of gas to our 
shock model results. Its corresponding emission was calculated 
with our LVG module used in a 'homogeneous slab' mode, and 
indeed generated significant emission only from the (3-2) and 
(4-3) lines, as can also be seen in Fig.|3] We used a linewidth of 
10 km s"' and a density of «h = lO'* cm"-' with a fractional abun- 
dance of X(CO) - 10"^ for the blue and red cases. We found 
that this 'CO layer' is warm (35 K and 15 K for the blue and 
red components, respectively), and adds to the pure-shock CO 
column density in both components, from 2x10'^ cm"^ for the 
blue component, and 10'^ cm"^ for the red one. Our final result, 
the sum of the C-shock and 'CO layer', is also shown in Fig. [3] 



To understand the physical origin of this emitting layer, we 
varied the cosmic ray ionisation rate in our best-fitting models 
from its solar value ^1 = 5 x 10"'^ s"' to the ^2 - 3.4 x 10"'^ s"' 
indicated bv iHewitt et alJ (l2009i) in the specific case of W28. 
In our modelling, this modification impacts on the chemistry 
(higher ionisation fraction) and the physics (warmer gas) of 
our outputs. The resulting shock profile is more naiTow (hence 
generating slightly less CO emission, see Fig. [3]l, with higher 
post-shock temperatures close to the aforementioned CO layer. 
Because one of the limitations of our shock model is the ab- 
sence of significant velocity gradient in the post-shock region, 
the warm CO layer component might correspond to a post-shock 
layer of gas mildly heated by the energetic radiation that seems 
to exist in the region (UV or more energetic, as characterised 
by iHewitt et al.l l2009r) . The rigorous inclusion of these radia- 
tive eiTects in our shock model is still in progress. In Sect. |B] 
we present a consistency check based on the use of the Spitzer 
H2 observations. 
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Fig. A.l. Location of the field covered by our CO observations 
on t he larger-scale radio c ontinuum image at 327 MHz taken 
from lClaussen et all (1 19971) : entire SNR in the left panel, zoom 
in the right panel. The shown CO observations are the (6-5) map, 
also displayed in Fig.[T] 

Appendix A: The APEX observations 

APEX observations towards the supernova remnant W28F were 
conducted in several runs in the year 2009 (in May, June, 
August and October). We used of a great part of the suite 
of heterodyne receivers available for this facility: APEX- 
2 Hisacher et al (2006), FLAS H 460 H e vmincketal (2006), 
and CHAMP+ iKasemann et aP (l20()6h : iGusten et alJ (llOOSh . 
in combination with the MPIfR fast Fourie r transform spec- 
trometer backend (FFTS, iKlein et aP l2006l) . The central po- 
sition of all observations was set to be q'[j2()oo] = 18''01"'5P.78, 
y6[j2000]=-23°18'58'.'50. Focus was checked at the beginning of 
each observing session, after sunrise and/or sunset on Mars, or 
on Jupiter. Line and continuum pointing was locally checked on 
RAFGL1922, G10.47B1, NGC6334-I or SgrB2(N). The point- 
ing accuracy was found to be of the order of 5" rms, regard- 
less of the receiver that was used. Table [ATI contains the main 
characteristics of the observed lines and corresponding observ- 
ing set-ups: frequency, beam size, sampling, used receiver, ob- 
serving days, forward and beam efficiency, system temperature, 
spectral resolution, and finally the velocity interval that was used 
to generate the integrated intensity maps. The observations were 
perform ed in position- s witch ing/raster mode using the APECS 
software iMuders et al.l ('2006). The data were reduced with the 
CLASS software (see http://www.iram.fr/IRAMFR/GILDAS). 
For all observations, the maximum number of channels avail- 
able in the backend was used (8192), except for CO (4-3), for 
which only 2048 channels were used, leading to the spectral res- 
olutions indicated in Table I A. fl Maps were obtained for all con- 
sidered transitions, covering the field introduced in Fig. [Tjand 
put in the perspective of the whole SNR in Fig. lA.ll 

Appendix B: The H2 observations 

B.1. The dataset 

As a consistency check for our models, we used of the 
SpitzerflRS observations of the H2 pure rotational transitions 
(0-0 S(0) up to S(7)), reported and analysed in N07 and Yll. 
Although the original dataset also includes other ionised species, 
we chose to use only the H2 data. The raw product commu- 
nicated to us by David Neufeld contains rotational transitions 
maps, with 1.2" per pixel, centred on q'[j2ooo] = 18''01 "'52^32, 
/5[j2000j=-23°19'24'.'92. Fig.|BlT]shows an overlay of our APEX 
CO (6-5) map with the H2 0-0 S(5) region observed by Spitzer. 
The figure shows coinciding maxima between the two datasets 
in the selected position, and a slightly different emission distri- 




offset (arcsec) 

Fig. B.l. Overlay of the map of CO (6-5) emission observed by 
the APEX telescope (colour background) with the H2 0-0 S(5) 
emission (white contours), observed with the Spitzer telescope. 
The wedge unit is K km s"' (antenna temperature) and refers 
to the CO observations. The H2 0-0 S(5) contours are from 50 
to 210 cr, in steps of 20cr - 1.6 x 10""^ erg cm"^ s"' sr"'. The 
green contour defines the half-maximum contour of this transi- 
tion. Like in Fig.[Tl the blue circle indicates the beam size of the 
SOFIA/GREAT observations, on the position of the centre of the 
circle, marked by a blue dot. The beam and pixel sizes of the CO 
(6-5), and H2 0-0 S(5) observations are also provided in yellow 
(big and small, respectively) circles in the lower right corner. The 
black contour delineates the half-maximum contour of the H2 0- 
S(2) transition. The field is smaller than in Fig.[T] and the (0,0) 
position is that of the Spitzer observations (see Appendix iBt. 



bution. This might be the effect of the better spatial resolution of 
the H2 data, which reveal more peaks than in CO. This overlay 
also shows the slightly different morphology of the emission of 
the S(2) (half maximum contour in black) and S(5) (half maxi- 
mum contour in green) transitions, at the available resolution. 



B.2. Excitation diagram 

We performed a consistency check on our modelling (see 
Sect. O using the excitation diagram derived for the selected 
emission region. The H2 excitation diagram displays ln{Nuj/gj) 
as a function of Euj/k^, where N^j (cm^^) is the column density 
of the rovibrational level (v, J), E^j/Icb is its excitation energy 
(in K), and gj - (2/ + 1)(2/ -H 1) its statistical weight (with 
7=1 and 7 = in the respective cases of ortho- and para- 
H2). If the gas is thermalised at a single temperature, all points 
in the diagram fall on a straight line. The selected position for 
the present study was introduced in Sect. 12.1] and can be seen 
in Figs. [U and IB. II The column density of the higher level of 
each considered transition was extracted by averaging the line 
intensities in a 1 1 .85" radius circular region, consistent with our 
handling of the CO data. In the process, we corrected the line 
intensities for extinction, adopting the visual extinction values 
fro m N07, A„ = 3 - 4, and u sing the interstellar extinction law 
of Rieke & Lebofskv (1985). The resulting excitation diagram 
can be seen in Fig. IB. 21 We note that the intensity of the H2 0-0 
S(6) transition (J = 8) cannot be determined reliably, because 
the line is blended with a strong 6.2 pm PAH feature. 



A. Gusdorf et al.: Probing MHD Shocks with high-7 CO observations: W28F, Online Material p 2 
Table A.l. Observed lines and corresponding telescope parameters for the APEX observations of W28F. 



line 


CO (3-2) 


CO (4-3) 


CO (6-5) 


CO (7-6) 


''CO (3-2) 


y (GHz) 
FWHM (") 
sampling (") 


345.796 

18.1 

10 


461.041 

13.5 

7 


691.473 
9.0 

4 


806.652 

7.7 
4 


330.588 

18.9 

10 


receiver 
observing days 


HET345 
04-05/08 


FLASH460 
04/06 


CHAMP+ 

13-14-15/08 


CHAMP+ 

13-14-15/08 


HET345 
05-06-07/08 


Feff 
Seff 


0.97 
0.73 


0.95 
0.60 


0.95 
0.52 


0.95 
0.49 


0.97 
0.73 


T,y, (K) 

Ao(kms"') 


279-288 
0.106 


372-439 
0.318 


1330-1722 
0.079 


3524-5850 
0.068 


357-377 
0.166 


reference offset (") 


(-100,-50) 


(-120,0)" 


(-100,-50) 


(-100,-50) 


(-100,-50)" 


u interval (km s"') 


-30/40 


-30/40 


-30/40 


-10/25 


-30/40 



Notes. *"' The observations were carried out with this off position, and were corrected a posteriori from its potential contamination, by means of 
an off spectrum on the (-120,0) position taken with an off position at (-240,0). 

'*' The observations were carried out with this off position, and were corrected a posteriori from its potential contamination, by means of an off 
spectrum on the (-100,-50) position taken with an oif position at (1000,1000). 



B.3. Comparisons witti our models 

Unlike our CO observations, the Spitzer/IRS observations are 
not spectrally resolved (owing to the relatively low resolving 
power of the short low module in the range 60-130). On the 
other hand, given the minimum energy of the levels excited in 
the H2 transitions (for 0-0 S(0), £„ ~ 509.9 K), we can make 
the double assumption that the lines are not contaminated by 
ambient emission or self-absorption, and that the measured line 
intensities are the result from the blue- and red-shock present in 
the line of sight. Therefore, one must compare the observed H2 
level populations to what is generated by the sum of our best- 
fitting blue- and red- shock models. Fig. lB.2l shows the compari- 
son between the excitation diagram derived from the sum of our 
two CO best-fitting models and the observed one. The H2 excita- 
tion diagram comprising the sum of our CO best-fitting models 
is shown in both panels in red circles, and provides a satisfying 
fit to the observations. With the aim of improving the quality of 
this fit, we also studied the influence of the initial ortho-to-para 
ratio (OPR) value, upper panel, and of the cosmic ray ionisation 
rate value, lower panel. 

In the upper panel, for a cosmic ray ionisation rate value 
adopted as the solar one, ^i = 5 x 10"'^ s"', the influence of 
the initial OPR value is studied. In their models, N07 inferred a 
value of 0.93 for the 'warm' component (~322 K). On the other 
hand, the value associated to their 'hot' (~1040 K) component 
could not be estimated, owing to the large uncertainty associ- 
ated to the determination of the H2 0-0 S(6) flux. In our shock 
models, the OPR value is consistently calculated at each point 
of the shocked layer, and is mostly the result of conversion re- 
actions between H2 and H, H^ or H^. The heating associated 
to the passage of the wave increases the OPR value towards the 
equilibrium value of 3.0. Nevertheless, this value is reached only 
when it is also the initial one. The N07 situation is then probably 
adequately approximated in our models where the OPR initial 
value is less than 3. However, the influence of this parameter is 
minimal on the excitation diagram, where it only seems to create 
a slight saw-tooth pattern between the odd and even values of J. 

We finally investigated the effects of the high-energy radia- 
tion field on the excitation diagram in the lower panel of Fig. IB. 21 
The limits of our modelling are reached because the only han- 
dle that we have to study this effect on our outputs is the vari- 
ation of the cosmic ray ionisation rate value, which affects the 
corresponding chemistry (higher ionisation fraction) and physics 



(warmer gas). A proper treatrnent of t he energetic radiation com- 
ponents listed in 'Hewitt et alj (l2009l) should indeed incorporate 
UV pumping or self-shielding, as would be the case in PDR re- 



gions (e.g.. lHabartetal.ll2011 



citation ef fects by X-ray (see 



), as well as chemica l and H2 ex- 
Dalgamo et al.lll999l) and cosmic 



ay I 

rays (e.g.. lFerland et al.ll2008l) . Keeping these limits in mind, we 
found that even strong modifications of the cosmic ray ionisa- 
tionrate from the solar value of 5x10"'^ s"' to the more extreme 
value of 10"'"^ s"' have only very limited effects on the excitation 
diagram. The rather convincing fits to the H2 data seem to indi- 
cate that these excitation effects by energetic photons could well 
be minimal in our case, but we repeat that their proper inclusion 
to our models is work in progress. 
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Fig.B.2. H2 excitation diagram comparisons. Upper panel: evo- 
lution of the modelled excitation diagram obtained for ([ - 
5 X 10 ''' s ', varying the initial value of the OPR from the un- 
realistic, extreme-case value to its equilibrium one, 3. Lower 
panel: influence of the cosmic ray ionisation rate variation on 
our modelled excitation diagram, with the initial OPR set to 3. 



